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ABSTRACT 



^^ ' Context. Active Galactic Nuclei (AGN) are thought to play an important role in galaxy evolution. It has been suggested that AGN- 

» I , feedback could be partly responsible for quenching star-formation in the hosts, leading to transition from the blue cloud to the red 

C^ . sequence. The transition seems to occur faster for the most massive galaxies, where traces of AGN activity have been already found 

at z < 0. 1 . The correlation between AGN activity, aging of the stellar populations and stellar mass still needs to be fully understood, 

especially at high redshifts. 

Aims. Our aim is to investigate the link between AGN activity, star- formation and stellar mass of the host galaxy at < z < 1, 
looking for spectroscopic traces of AGN and aging of the host. This work provides an extension of the existing studies at z < 0. 1 and 
contributes to shed light on galaxy evolution at intermediate redshifts. 

Methods. We used the zCOSMOS 20k data to create a sample of galaxies at z < 1 . We divided the sample in several mass-redshift 
bins to obtain stacked galaxy spectra with an improved S/N. We exploited emission-line diagnostic diagrams to separate AGN from 
star-forming galaxies. 

Results. We found indication of a role for the total galaxy stellar mass in leading galaxy classification. Stacked spectra show AGN 
signatures above the log M,/Mq > 10.2 threshold. Moreover, the stellar populations of AGN hosts are found to be older with respect 
to star-forming and composites galaxies. This could be due to the the tendency of AGN to reside in massive hosts. 
Conclusions. The dependence of the AGN classification on the stellar mass is in agreement with what has been already found in 
previous researches. It is consistent with, together with the evidence of older stellar populations inhabiting the AGN-like galaxies, 
the downsizing scenario. In particular, our evidence point to an evolutionary scenario where the AGN-feedback is capable of quench- 
ing the star formation in the most massive galaxies. Therefore, the AGN-feedback is the best candidate for initiating the passive 
evolutionary phase of galaxies. 

Key words. Galaxies: active - Galaxies: evolution - Galaxies: starburst - Galaxies: stellar content 



1. Introduction 

In recent years, a better overview of galaxy formation and 
evolution has been developing, with a great contribution coming 
from large cosmic surveys. In particular, the local galaxies 
colour function has b een discovered to be a bimodal function 
dStrateva et all l200lt 



I2004t iBalogh et alj|2004l: iHogg et al.ll20(> 



KauffmannetaP l2003t 



iBaldrv et all 



The two dominant 
sequences are associated with the ellipticals+SOs population 
(red distribution, or red sequence) and the spirals+irregulars 
population (blue distribution, or blue cloud), while galaxies 
showing intermediate properties are considered as composites or 
transitional objects (green valley). Early-type galaxies (ETGs) 
dominate at the top-end of the local galaxy mass function 
where, on the contrary, late-type galaxies (LTGs) become 
progressively more rare. The discovery of how the black hole 
(BH) ma ss strongly correlates with the properties of the ho sting 
galaxies (iMagorrian et aI.III998T: iFerrarese & MerrittlEoOOl ). has 
generated new hypotheses in linking the black hole growth with 
galaxy formation and evolution. The mass of the host seems 
to play a key role in the lat ter process as well . According to 
the " d ownsizing" scenario ( Cowie et al. Il996t iKodama et al. 
20041: iThomasetalJ 120051: JTanaka et alj 120051: iJuneau et al. 
20051: iBundv et alj l2006h . the sites of active star-formation 
include high-mass galaxies at high redshift (z > 1) but only 
intermediate- and low-mass galaxies at z < 1. This scenario 
is supported by several studies, e.g. it has been found that 
luminous and massive spheroidals build up most of their stellar 
mass prior to z ^ 1, while l ow mass ellip t icals c ontinue to 
grow after this cosmic time (iHeavens et al.l 20041: iTreu et all 
l2005t Ivan der Wei et al.ll2005L iJimenez et al.ll2005h . Downsizing 
may be a natural expectation in a hierarchical galaxy formation 
scenario, provided that there is a mechanism that quenches star 
formation at earlier times in massive galaxies. As a support to 
downsizin g and the exi s tence of a link between mass and galaxy 
evolution. iBundv et alj ( |2004|) found that the relative abundance 
of red galaxies with intermediate mass (M a; 6 x 10 "'M©) 
increases by a factor of 3 from z ~ 1.2 to z ~ 0.55, while the 
number of blue late-type galaxies declines significantly. 
The fast suppression of the star formation activity (quenching) 
is currently thought to initiate passive galaxy evolution and to 
be largely respon sible for the growi ng abundance of galaxies on 
the red sequence dFaber et al.ll2007l) . The quenching mass (Mq), 
namely the stellar mass above which galaxies start their passive 
phase , has found to be a function of the redshift dBundv et al.l 
l2004l) . Blue and star-forming galaxies at their quenching mass 
are redder than the average blue population, suggesting a 
transion from blue- to red-type galaxies. 

The star-formation quenching could be explained by several 
mechanisms. The observed decline of the galaxy star formation 
from z ~ 1 to z - can be due to a gradual gas exhaustion 
(ICooper et al.l l2008h . Another mechanism that could explain 
the quenching is major mergers. The interaction between a 
galaxy and its neighbours could indeed produce gravitational 
torques on the gas, reducing its angular momentum and sending 
it toward the gala ctic center after an initial starburst phase 
(JMartig et al.ll2009l) . As a result, fuel is available for accretion, 
allowing the growth of the central Black Hole (BH) rather than 
the formation of new stars in the disk. However, it has been 
established that mergers cannot be the main cause to the build 
up of the galaxy bimodality, since luminous and massive old 
galaxies formed throughout this process were already common 
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at z ~ 1 (IConselice et al.ll2007l) and their number density only 
declines above this value. If merger events were efficient in 
forming galaxies at z < 1, we would expect an increasing 
number of massive objects at low redshift, which has not been 
observed. Similarly, a role for active galactic nuclei (AGN) 
has been suggested (ICiotti & Ostrikej Il997t iBinnevI 12004 : 
Isii3l2005t ISpringel et al.ll2005l) . BHs are thought to be a basic 
constituent of most massive systems (JRichstone et al. 1998]) 
and they have also been identified in some late-type and dwarf 
galaxies (JFilippenko & Hoi l2003l) . During the BH accretion, 
AGN release a large amount of energy, which could be partially 
absorbed by the host galaxy to quench the star formation 
and eventua lly transform blue g alaxie s into red galaxies (e.g . 
ISilk&Reeslll998l: ICroton et al.l 120061 iNaravanan et al.l llool . 
According to some calculations (ICattaneo et al.l 120091) . even a 
small fraction (< 1%) of the energy released within a bulge 
during the BH accretion would be sufficient to heat and/or blow 
away the entir e gas c ontent. Consistently with this framework , 
iHopkins et aP (l2006l) further developed ISanders et all (Il988h 
scenario in which starburst, quasar active phase, BHs growth 
and elliptical galaxies are connected to each other in a galaxy 
evolutionary sequence and the gas consumption is coupled with 
supernova-driven winds and/or AGN-feedback to quench the 
star formation. 

The hypothesis of the AGN-feedback is strengthened by the fact 
that a large amount of galaxies have been found to go through 
an active phase. Summed over all Hubble ty pes, roughly hal f 
of all galaxies can be considered as active dHo et al.lll997al) . 
even though some of them show a fairly low activity level. 
The amount of active galaxies becomes more remarkable for 
galaxies with a prominent bulge component, rising to 50 - 70% 
for Hubble types E-Sb, whereas the detection rate of AGN 
drops towards later Hubble types (Sc and later), where about 
80% of the galaxies host a star-forming nucleus instead. This 
sugg ests a clear de pendence of nuclear activity on the Hubble 
type (iKauffmann et al. 20031 iMiller et al.ll2003l) . The observed 
distribution of Hubble types for galaxies hosting active versus 
inactive nuclei leads to the expectation that the blue and red 
populations could have not only different nuclear properties but 
global as well. In particular, it seems that the observed high 
fraction of AGN in the green valley points to the connection 
between the end of the star formation phase, observed in a 
change of the galaxies color, and the rise of the AGN activity 
(Silverman et al. 2009; Schawinski et al. 2010). 
AGN can be identified by a variety of methods: an unusually 
blue continuum, strong radio or X-ray emission, the presence of 
strong or broad emission-lines. In nearby galaxies hosting low- 
luminosity AGN, we expect the non-stellar signal of the nucleus 
to be weak with respect to the one coming from the host galaxy, 
thus one of the best and least biased methods to search for AGN 
is to conduct a spectroscopic survey of a complete, optical-flux 
limited sample of galaxies and study their spectral properties. 
The line-emitting gas in AGN and star-forming galaxies (SFGs) 
is powered by two different ionizing mechanisms, producing 
different emission line ratios: accretion around black hole and 
photoionization by hot massive OB stars. It is important to note 
that emission lines in narrow-line AGN have a considerably 
great ionization range and that low-ionization lines are stronger 
than in normal star-forming galaxies. It has been demonstrated 
that ordinary O-type stars do not produce sufficiently strong 
low-ionization lines, accounting for the observed spectra, as 
they onl y account for a fe w percentage of the total blue light 
at most (ISarzi et al.ll2005h . Therefore, photoionization arising 
from powerful central non-thermal source is the ideal candidate 



M. Vitale et al.: AGN in the zCOSMOS survey 



to explain the excitation mechanism in galaxi es which display 

strong low-ionization emission-lines (lHdl2005l) . 

Emission-line dia gnostic diagr a ms ( Baldwin et al.l 



2. Data 



198U IVeilleux & Osterbr oca Il987l iKewlevet alJ l2003t 



LamareilieetalJ 12004; .Lamareillell2010t) represent a powerful 2. t. COSMOS and zCOSMOS surveys 



tool for probing the nature of the dominant ionizing source in 
galaxies, hence to distinguish between objects dominated by 
star formation and galaxies where the nuclear activity is more 
relevant. Narrow-line AGN can be identified in spectroscopic 
galaxy surveys by the ratio of some distinctive emission-lines, 
such as [NII]^6583 A (hereafter [Nil]) or [SII]^^6717,31 A 
(hereafter [SII], referred to the combined luminosity of the 
doublet) over HQ;i6563 A and [OIII] A5001 A over H^/14861 
A (hereafter Ha, [OIII] and HjS, respectively). The lines that 
appear in the reddest part of the spectrum are shifted out of the 
visible wavelength range for z > 0.5. In this case, it is possible 



to use the 



Oll]AA3726 + 3729 A l i ne (hereafter [Oil]) instead 



('Rolaetal.f ll997t iLamareille et al.l l2004t iPerez-Montero et al.l 
[2007), which is optically visible up to z ~ 1. 
The aim of this work is to investigate whether it is possible 
to highlight some evidence of AGN activity in suppressing 
star-formation, with a particular look on a possible trend of 
the galaxy total stellar mass. This attempt has been already 
made in the past at low-redshift (z < 0.1) and the analysis 
is now extended to redshift up to z ~ 1. iKauffmann et alj 
(I2003h find, by using Sloan Digital Sky Survey (SDSS) data, 
that galaxies divide in two distinct families at a stellar mass 
threshold of 3 x 10'" Mq. The least massive galaxies show 
young stellar populations and the low concentrations which 
are typical of LTGs. As the stellar mass increases, galaxies 
show older stellar populations and higher concentrations, 
typical of ETGs. By using emission-line diagnostic diagrams, 
they show a dependence of the AGN-detection rate on mass, 
although the authors warn that many objects are classified as 
"composites" (contributions to the emission-line spectrum come 
by both star-formation and AGN) because of the large size of 
the SDSS fiber, which can collect up to 40% of the total light 
from the galaxy. Blue star-forming objects have been found 
to have masses at the low end of the mass function, while the 
majority of the ETGs lie at the high-end of t he mass function 
( jBaldrv et al.1 1200 4'; 'Schawins ki et al.l |2007|) . Very recently, 
however. JAird et a l. (2012) argued that the finding of more AGN 
in massive hosts is due to selection eff'ects. The latter seem to be 
driven by the Eddington-ratio distribution of AGN in galaxies 
of any stellar mass. In particular, AGN are more easily detected 
in massive hosts as sources accreting at low Eddigton rate are 
more luminous than in less massive hosts. 
We analyze the evolution of galaxy spectral properties up to 
z ~ 1 by using ga laxies from the zCOSMOS-Bright 20k sample 
(iLillv et al.l2007h . With the large number(~ 2x10^*) of observed 
objects, we are able to measure emission lines from high S/N 
stacked spectra. The paper is organized as follows. In Section 
2 we present the zCOSMOS survey and the sample selection. 
In Section 3 we present the spectral analysis of the stacked, 
stellar-continuum subtracted galaxy spectra. In Section 4 we 
exploit a set of classical and more recent diagnostic diagrams 
to identify the main ionizing process inside galaxies and look 
for a trend of the stellar mass. Section 5 contains our results 
and discussion. Main findings and conclusions are presented in 
Section 6. 



The COSMOS survey is a large HST-ACS surve y, with /-band 
expos ures down to Iab = 28 on a field of 2 deg^ dScoville et al.l 
120071) . The COSMOS field has been the object of extensive 
multiwavelength ground- and space-based observations span- 
ning the entire spectrum: X-ray, UV, optical/NIR, mid-infrared, 
mm/su bmillimeter and radio , providing fluxes measured over 30 
bands (iHasinger et al.ll2007l: Taniguchi et al. 20071: Caoak et al.l 



'2007/Lillv et al."2007'; 'Sanders et al."2007';'Bertoldi et alj|2007t 
Schin nerer et al. 2007; Koekemoer et al. 2007; McCracke n et al.l 
120101) . 

The zCOSMOS survey was planned to provide high-qua lity red- 
shift information to the COSMOS field ilLiUv et alj2"007l) . It ben- 
efitted of ~ 600 h of observations at VLT using the VIMOS 
spectrograph and it consists of two parts: zCOSMOS-bright and 
zCOSMOS-deep. The zCOSMOS-deep targets ~ 10 000 galax- 
ies within the central 1 deg^ of the COSMOS field, selected 
through color criteria to have 1.4 < z < 3.0. The zCOSMOS- 
bright is purely magnitude-limited and covers the whole area of 
1 .7 deg2 of the COSMOS field. It provides redshifts for ~ 20 000 
galaxies down to Iab S 22.5 as measured from the HST-ACS 
imaging. The success rate in redshift measurements is very high, 
95% in the redshift ran ge 0.5 < z < 0.8 , and the velocity accu- 
racy is ~ 100 km s ' jLillv et al.ll2009l) . Each observed object 
has been assigned a flag according to the reliability of its mea- 
sured redshift. Classes 3.x, 4.x redshifts, plus Classes 1.5, 2.4, 
2.5, 9.3, and 9.5 are considered a secure set, with an overall reli- 
abihty of 99% (see Lilly et al. 2009, for details). 
Our work is based on the the zCOSMOS-bright survey final re- 
lease: the so called 2Qk sample, totaling 16 623 galaxies with 
z < 2 and secure redshifts according to the above flag classifica- 
tion. It includes 18 206 objects in total, including stars. 
For objects brighter than Iab = 22.5 and without secure spec- 
troscopic redshift, photometric data from the COSMOS surve y 
provides good quality photometric redshifts (lllbert et al.ll2009l) . 
Based on a compa r ison w ith the zCOSMOS spectroscopic red- 
shifts, lllbert et aP (l2009l) estimate an accuracy of cr~pi,„, - 
0.007 X (1 -H Zj) for galaxies brighter than Iab =22.5. 
For all galaxies brighter than Iab = 22.5, absolute rest-frame 
magnitudes and stellar masses were obtained using standard 
multi-color spectral energy distribution (SED) fitting techniques, 
using the secure spectroscopic redshift, if available, or the photo- 
metric one^_Stell^imsseswere_obtah^ hyperzmass 
code dPozzetti et al.ll2010t iBolzonella et al. l2010l) . by assuming 
a Chabrier initial mass function (IChabrierll2003b . To establish 
mass completeness. iPozzetti et all (1201 Ol) define at each redshift 
a minimum mass (M„„„) above which the derived galaxy stellar 
mass function is essentially complete (all types of galaxies are 
potentially observable). To derive M^in, the limiting stellar mass 
(Miiiii) of each galaxy is calculated, namely the mass the galaxy 
would have, at its spectroscopic redshift, if its apparent magni- 
tude were equal to the limiting magnitude of the survey. Then, in 
order to derive a representative limit for the sample, the authors 
use the M/,„, of the 20% faintest galaxies at each redshift. M,„i„{z) 
is defined as the upper limit of the M/„„ distribution below which 
95% of the M/„„ values at each redshift are represented. This is 
assumed to be the completeness limit of the galaxy stellar mass 
function. 
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2.2. Analyzed galaxy sample 

In order to create the bins that collect the galaxies to be stacked, 
we used the Bright zCOSMOS spectroscopic catalogue {v4.12 
version). In case of AGN that are more luminous than the host 
galaxy, as the type-1 AGN where the emissio n is not screened 
by the dusty torus (AGN-unification model by iPadovani & Urrvl 
II992I) . we are not able to see the emission lines from the host 
galaxy. Therefore, we remove type-1 AGN to avoid overshin- 
ing issues. This was performed by excluding all the objects with 
broad-emission lines. Furthermore, we decided to exclude from 
our sample all the sources that do not bring an indication of 
well measured spectroscopic redshift. These latter classes and 
the broad-line AGN constitute, all together, around 6% of the 
original sample within 0.1 < z < 1 (16 678 objects). In particu- 
lar, we exclude 892 low-flag objects and 35 broad-line AGN. The 
galaxies with reliable flag are then 15 715. Another significant 
fraction of galaxies is lost during the binning process due to ex- 
clusion of underrepresented regions of the redshift-stellar mass 
diagram (see next section). After applying these selection crite- 
ria, we are able to study roughly half of the original zCOSMOS - 
Bright sample. 



3. Spectral analysis 

In this chapter we present how we exploited the large number 
of zCOSMOS galaxies by dividing them in various redshift and 
mass bins. This was done in order to obtain average spectra with 
high S/N and to investigate the evolution of the spectral proper- 
ties of the stacked galaxies. We aim to study the variation of av- 
erage properties of galaxies that are representative of a particular 
mass-redshift bin. After obtaining the spectra, we estimate and 
subtract the stellar contribution. Then, we measure line fluxes 
and compute the visual extinction. Finally, we use diagnostic 
diagrams to distinguish between the main ionizing mechanism 
inside these objects. Several programs and tasks have been used 
for this work, most of them running in the Image Reduction and 
Analysis Facility (IRAF) environment. 



3.1. Sample binning 

In order to search for variations of galaxy spectral properties as 
a function of redshift and stellar mass, we divided our large sam- 
ple in bins. The bin width was tested in order to achieve a good 
signal to noise with the largest number of mass-redshift bins. In 
particular, we tested the increase in S/N according to the increase 
in the number of stacked spectra. The first indication of a link be- 
tween total galaxy stellar mass and classification in the spectral 
types - AGN, composite galaxies, star-forming galaxies - made 
us decide for a fine grid in M» . On the contrary, the variations we 
applied to the size of the redshift intervals did not affect the clas- 
sification significantly. We afterwards decided for the number of 
z bins to be equal to 8, with a redshift step ~ 0. 1 dex. We divided 
the sample in 9 mass bins, with a mass step of 0.25 dex. The 
bins include few hundreds of objects for each redshift and mass 
range (See Tab.lTJ. For our sample selection, mass completeness 
is achieved when including in the bin s galaxies above the com- 
pleteness curves (iPozzetti et al.ll2010h . While at lower redshift 
most bins lie above the completeness limits, at higher redshifts 
roughly half of the bins are placed below the confidence curves. 
Although, in this case, the sample suffers from incompleteness, 
we chose to study all the high-redshift bins in order to probe a 
comparable dynamical range of masses to the low-redshift case. 



The choice of the bins comes from the following further consid- 
erations: 

- below 9 and over 11.25 logM^/M©, objects are poorly rep- 
resented in our sample. The S/N ratio of the average spectra 
made out of the stacking of an insufficient number of galax- 
ies would not allow precise line measurements; 

- 8 bins (4 per region covered by each of the two types of di- 
agnostic diagrams, z < 0.5 for the low-redshift diagnostic 
diagrams and z > 0.5 for the high-redshift diagnostic dia- 
gram) are a good trade off between the need to explore the 
cosmological evolution of the galaxy spectral properties and 
to keep the number of objects per bin (thus the S/N) high 
enough. 

For diagnostic diagrams based on [Nil], [SII], Ha, [OIII] and 
Hj6, the affected rest-frame spectral range is from 4800 to 6650 
A. The zCOSMOS spectra cover the 5650 to 9550 A range, 
after cutting off 100 A at both extremes to avoid noisy re- 
gions where the flux calibration is more uncertain. The useful 
redshift range for the low-redshift diagnostic diagrams is then 
0.177 < z < 0.436. In the same way, the redshift range for the 
high-redshift diagnostic diagram is: 0.548 < z < 0.884. 
In Fig. [T] we plot the zCOSMOS galaxies as a function of the 
total galaxy stellar mass. The plot displays our bin selection 
and the large scale structure within the zCOSMOS field. The 
general trend is a clear increase of the average stellar mass 
with the redshift due to the well-known selection effects of a 
magnitude-limited survey, with overdense regions (filaments), 
e.g. at 0.3 > z > 0.4 or 0.7 > z > 0.8, and underdense regions 
(voids), e.g. at z ~ 0.25. At low redshift the minimum mass 
{M„,i„) of the global galaxy population is better represented by 
the blue population, whereas it is closer to the red population if 
we are at higher redshift. This is because blue galaxies dominate 
the mass function at low redshift, due to selection effects. On the 
other hand, at high redshift the global M„„„ shifts towards the 
limit of the reddest population, since the galaxies have a higher 
mass-to-light ratio (M/L) on average. 

3.2. Stacking 

After the binning process, the spectra were stacked. Since each 
stacked spectra is the result of the average of up to several hun- 
dreds of single galaxy spectra included in a single bin, the gain 
in S/N is considerable and allows a more accurate line fitting. In 
order to create the composites, each spectrum was shifted to the 
rest-frame. The spectra were normalized in wavelength ranges 
always present in the observed spectroscopic window and lack- 
ing prominent spectral features. The wavelengths of the chosen 
ranges are: 5300-5800 A for 0.177 < z < 0.371; 4500 -4800 A 
forO.371 < z < 0.716; and 3400-3700 A for 0.7 16 < z < 0.884. 
Stacked spectra are shown in Appendix A. The main emission- 
and absorption-lines are labelled on top of each figure. 
The stacked zCOSMOS spectra present, especially for the ob- 
jects with low masses, a quite blue and featureless continuum 
at low redshift. The spectra are characterized by the presence 
of strong low-ionization emission lines such as [Nil], [SII] and 
[OIII] and strong recombination lines of the Balmer series, as 
expected from the spectral contribution of a young stellar pop- 
ulation. On the other hand, the continuum is much redder and 
displays a strong D4000 break at higher redshift, which is in- 
dicative of an older stellar population. The spectra show that the 
most massive galaxies are spectroscopically more evolved at at 
a given redshift, according to the downsizing scenario. 
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Fig. 1. The stellar mass of zCOSMOS 20k galaxies selected for our study (secure flags, no broad-line AGN) plotted as a function 
of redshift. The blue dashed lines represent, respectively, the mass completeness limits of the global galaxy population (solid), 
the population of early-type (dashed) and the late-type galaxies (dot-dashed). The black boxes mark the 61 mass-redshift bins that 
define the stacked galaxies we study with the diagnostic diagrams. 



Table 1. Number of galaxies in each mass-redshift bin of the sample. 



M„l Zm 


0.21 


0.27 


0.34 


0.4 


0.59 


0.67 


0.76 


0.84 


11.12 


- 


24 


83 


47 


52 


133 


101 


114 


10.87 


51 


42 


148 


107 


106 


241 


237 


189 


10.62 


63 


45 


196 


118 


146 


324 


296 


221 


10.37 


68 


69 


177 


143 


168 


323 


272 


198 


10.12 


57 


66 


189 


156 


187 


283 


216 


173 


9.87 


52 


79 


214 


152 


179 


286 


219 


201 


9.62 


76 


73 


213 


162 


193 


238 


185 


87 


9.37 


83 


95 


266 


211 


- 


- 


- 


- 


9.12 


108 


109 


- 


- 


- 


- 


- 


- 



Notes. The table reflects the order of the mass-redshift bins shown in Fig. 1. The median mass (first column) is expressed in log M,/Mq. 
median redshift is indicated on top of each row. 



The 



3.3. Stellar continuum subtraction 

Galaxy emission lines - in particular the ones from the Balmer 
series - are often contaminated by absorption lines produced by 
the atmospheres of fairly massive stars within the galaxies. Since 
our spectra are characterized by both emission from the hot ion- 
ized gas and absorption from the stellar component, we need 
an effective strategy to make precise measurements of spectral 
features involved in the diagnostic diagrams. The technique we 
adopted to remove the starlight from an integrated spectrum is 



to subtract a suitable model of stellar continuum to each stacked 
galaxy spectrum to obtain an almost pure emission-line spec- 
trum. This is particularly helpful in correcting the strength of 
Balmer lines, which are involved in both the AGN identifica- 
tion and the computation of the visual extinction. The proce- 
dure also aims to additional information on stellar populations, 
like age and metallicity, that would have been lost if we had 
focused on emission lines only. For this purpose, we have devel- 
oped an user- written IRAF task that performs the spectra contin- 
uum subtraction. The spectral library to be fitted to the stacked 
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Fig. 2. Four examples of stellar continuum subtraction at low redshift (z < 0.436, one per redshift bin). For each panel, the original 
stacked spectrum (green), the most suitable stellar continuum to be subtracted (orange) and the resulting zCOSMOS spectrum after 
the stellar continuum subtraction (black) are shown. The main emission and absorpion lines are labelled on top. From top to bottom 
and from left to right both total stellar mass and redshift increase. Top-left panel: stacked spectrum of the bin 9.25 < log Mt/M^ < 
9.5, 0.177 < z < 0.242, stellar template representing a 290 Myr old simple stellar population (SSP) with Z=0.2 and E(B-V)=0.4; 
top-right panel: 9.75 < log M^/Mq < 10, 0.242 < z < 0.306, 12 Gyr old stellar template with exponential decay of the star 
formation, Z=0.2 and E(B-V)=0.2; bottom-left panel: 10.25 < log M./Mg < 10.5,0.306 < z < 0.371, 900 Myr old SSP with Z=0.5 
and E(B-V)=0.4; bottom-right panel: 10.75 < log MJMq < 11, 0.371 < z < 0.436, 11 Gyr old SSP with Z=0.2 and E(B-V)=0.4. 
The stacked spectra before stellar continuum subtraction and the stellar templates have been shifted up of arbitrary quantities to 
offer a clearer view. 



spectra must contain enough information on various absorption 
features, in order to be able to simulate the stellar components 
of galaxies within a wide range of ages and metallicities. We 
adopted the Bruzual and Chariot (BC) library, which is a library 
of stellar population synthesis models computed by Bruzual and 
Chari ot using their Isochrone S ynthesis Spectral Evolutionary 
Code (iBruzual & Charlotll2003l hereafter BC03). This code al- 
lows to predict the spectral evolution of stellar populations in 
various ranges of ages and metallicities at a resolution of 3 A, 
across the whole wavelength range (from 3200 A to 9500 A). 
Models with different time scales of star formation, initial mass 
functions and metallicities well reproduce the spectral and pho- 
tometric properties of nearby galaxies with various morpholog- 
ical types, from young irregulars to elliptical galaxies. For each 
BC template, we created 5 spectra with increasing B(E-V) val- 
ues (0.2, 0.4, 0.6, 0.8 and 1) in order to take reddening due to 
dust into account. A total of 234 stellar templates were com- 
pared to the stacked spectra to find the most suitable to be sub- 
tracted to each composite. The continuum subtraction task de- 
grades the resolution of the template spectra to the one of zCOS- 
MOS (R~ 600). Then, the template is chosen to have the low- 
est residuals in pre-selected regions of the continuum subtracted 
sample which are free from strong emission lines. These spectral 
windows have a typical width of ~ 100 A. They span the range 
between 3600 and 7300 A, for an overall coverage of about 2300 
A. 
The stellar continuum subtraction procedure successfully cor- 



rects for the underlying stellar absorption (Fig. 13). This fea- 
ture is particularly noticeable in the highest redshift bins (Fig. 
[3]l, where the strong stellar contribution to the galaxies spectra 
comes from older stellar populations. 



3.4. Emission line measurement 



We measured the flux of [Oil], Hy^4340 (hereafter Hy), H/3, 
[OIII], Ho-, [Nil] and [SII]. We used the Splot interactive IRAF 
task to make spectral measurements of gaussian-like emission 
lines. Ha and H/3, along with the oxygen lines, [Nil] and [SII], 
are used in the exploitation of the diagnostic diagrams. The er- 
rors on the measurements of the emission lines in our stacked 
spectra corrected for the stellar absorption are typically less 
than 20%, estimated on the basis of repeated independent mea- 
surements. The errors mainly depend on the line intensity and 
whether the emission line is blended with other spectral features. 
The resolution of the zCOSMOS spectra is sufficient to resolve 
the [SII] doublet in most of the composites, but it is not sufficient 
to resolve the [Oil] doublet, which always appears as a blended 
feature. 

For the highest stellar masses, the emission line measurements 
are more difficult because the S/N ratio is lower due to poorer 
statistics, besides the fact that the lines become weaker and the 
continuum stronger, especially at high redshift. 
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Fig. 3. As in Fig. IH but at higher redshift (z > 0.548). Top-left panel: stacked spectram of the bin 10.5 < log MJMq < 10.75, 
0.548 < z < 0.632, stellar template representing a 1.4 Gyr old SSP with Z=0.08 and E(B-V)=0.2; top-right panel: 10.75 < 
logM,/Mo < 11, 0.632 < z < 0.716, 1.4 Gyr old SSP with Z=0.2 and E(B-V)=0.2; bottom-left panel: 10.75 < log M./Mg < 11, 
0.716 < z < 0.8, 1.4 Gyr old SSP with Z=0.08 and E(B-V)=0.2; bottom-right panel: 11 < logM./Mo < 11.25,0.8 < z < 0.884, 1.4 
Gyr old SSP with Z=0.2 and E(B-V)=0.2. The stacked spectra before stellar continuum subtraction and the stellar templates have 
been shifted up of arbitrary quantities to offer a clearer view. 



3.5. Reddening correction 

Observed fluxes need to be corrected for extinction. The total 
visual extinction. Ay, can be determined by using the ratio 
between two recombination lines of the Balmer series. For this 
purpose. Ho- and HyS are commonly preferred as they are strong 
lines placed in an easily accessible spectroscopic position. In 
this work, we assumed th e Ha/H /J theoretical ratio to be 2.86 
and the Hy/Hd to be 0.47 (lOsterb rock 1989, case B recombina- 
tion, temperature T= 10'* K, electron density n^ - 10^ cm~^) 
and we assumed the Calzetti extinction curve (ICalzetti et al.l 
[T99l . with Rv = 4.05 + 0.8. 

In diagnostic diagrams that make use of lines which are close 
to each other in the spectrum, their ratio is relatively insensitive 
to reddening effects. Nevertheless, we chose to correct for the 
visual extinction because the [Oil] and Hfi lines are quite far 
away from each other in wavelength, thus their ratio is sensitive 
to reddening. 



The redshift range 0.436 < z < 0.548 has been excluded from 
our study because [Nil] and Ha are redshifted outside the optical 
spectra and [Oil] is not yet visible in the blue part of the spectral 
range. Therefore, it is not possible to use the first two diagnostic 
diagrams. In case of the MEx diagram, we decided to keep the 
same grid of stellar mass and redshift bins. Although it would be 
technically possible to include the missing redshift range when 
using the MEx diagram, we expect that this would not change 
our finding that the AGN classification is more dependent on the 
stellar mass rather than on the redshift. 



4. Diagnostic diagrams 

In this work, we adopted the following emission-line diagnostic 
diagrams: 



- [NII]/Hg versus [O III]/H/J (z < 



Il/Hor 

2J^ 



198U VeiUeux & OsterbrocM 1987 



c^ 



Kauffman n et al.ll2003tlKewlev et allboO^ 



- [OIIl/H/g versus [OIII]/H;g (z 



iRola et al Jl 1 9971: iLamareiUe et al 



- MEx diagram d Juneau et al.ll201 ll) 



4.1. Low redshiU diagnostic diagrams (Q.lll <z< 0.436J 

In the [NII]/Hg versus [OIII]/Hy8 or Baldwin-PhilHps-Terievich 
(BPT) diagram ([Baldwin et alJ[l98ll) . galaxies are distributed in 
two arms. The increase of the [NII]/!!^ ratio is a linear func- 
tion of the nebular metallicity and it presents a saturation point 
(Denicol, Terlevich & Terlevich 2002; Kewley & Dopita 2002, 
Pettini & Pagel 2004) above which any further increas e in the 
[NII]/Ha value is only due to AGN contribution (Kewl ev et al.l 
120061 : IStasihska et alJl2006r) . The BPT emission-line diagnostic 
diagram makes use of difl^erent demarcation curves, both theoret- 
ical and based on observations. The first curve was derived the- 
oretically by Kewley et al. (2001) in order to find an upper limit 
for star-forming galaxies (Fig.|4}. Kauffmann et al. (2003), using 
the large sample of emission line galaxies in the SDSS, defined 
.436) ([Baldwin et alJ a demarcation which traces more closely the observed lower left 
branch attributed to purely star forming galaxies. This results 

in a larger fraction of galaxies residing on the AGN side (Fig. 

0.548) (ITresse et aL 1996t [Hi. According to the latter classification, some of the AGN and 

composite galaxies with widely spread distributions of metallic- 
ity and ionization parameters are degenerate with star-forming 
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Fig. 4. Classical [Nil] diagnostic diagram. In both panels, a theoretical demarcation curve (dashed, Kewley et al. 2001), separates 
star-forming galaxies and composites from AGN, while an empirical demarcation curve (three-dot-dashed, Kauffmann et al. 2003) 
separates pure SFGs from composites and AGN. Each color represents the trend with stellar mass for a fixed redshift bin: 0.177 < 
z < 0.242 (blue), 0.242 < z < 0.306 (green), 0.306 < z < 0.371 (purple) and 0.371 < z < 0.436 (red). The total stellar mass 
increases from left to right. In the left panel, the SDSS galaxies from lVitale et al.l (1201 2h are plotted in gray. The right panel offers a 
a closer view of the zCOSMOS galaxies. 
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Fig. 5. [Oil] diagnostic diagram. [OII]/Hy6 fux dust-corrected ratios are used. The solid demarcation lines are bv iLamareillel (|20 1 Ol) . 
Each color represents the trend with stellar mass for a fixed redshift bin: 0.548 < z < 0.632 (blue), 0.632 < z < 0.716 (green), 
0.716 <z < 0.8 (purple) a nd 0.8 <z< 0.884 (red). The total stellar mass increases from left to right. In the left panel, the SDSS 
galaxies from lVitale et al.l (|20I2|) are plotted in gray. The right panel shows a closer view of the zCOSMOS galaxies. 



galaxies. This is especially the case, for example, of the [SII]/Ha 
versus [OIII]/Hj6 diagram. 

Objects that are classified as AGN in the [NII]-based diagram 
can be classified as star-forming galaxies in the [SII]/!!!! versus 
[OIII]/H;6 diagram. This is thought to be related to the enhance- 
ment of [SII] lines in starburst galaxies, due to the mechanical 
energy released into the gas by supernovae and stellar winds 



(also called shock excitation. iDopita et al.ll2002h . However, it is 
still not clear whether the [SII] enhancement is driven by highly 
ionizing photons produced by accretion onto the super massive 
black hole or if the SNe winds are a sufficiently powerful source 
of ionization. The [SII] -based diagnostic diagram, due to the 
fringing and the lines blending, is more affected by errors and 
leads to a higher probabiUty of misclassifying galaxies with re- 
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spect to the [NII]-based diagram. For these reasons, the diagnos- 
tic diagram that makes use of the [SII] doublet is not used in this 
paper to draw conclusions but it is only shown later as a compar- 
ison. 

In Fig.|4](left panel), we show the zCOSM OS data superimposed 
to the SDSS data from i Vitale et al.l(l2012l) . Each color represents 
a different redshift bin and each point along the curve repre- 
sents a different mass-bin. The right panel offers a zoom on the 
region of the diagram where our data are placed. The stacked 
spectra show a progressive shift towards the AGN region - on 
the right-hand side of the diagram - for increasing stellar mass 
(from left to right). The objects with the highest mass are placed 
on the right of the Kauffmann demarcation curve, which sep- 
arates the region where galaxies start to show the sign of sig- 
nificant nuclear activity. The stacks follow, for a fixed redshift 
bin, traks that move from the upper part o f the SF galaxies se - 
quence, down to the metallicity sequence dKewlev et al .1120061) . 
For log M,/Mq > 10.2, the objects start being classified as com- 
posites or AGN, in agr eement with what already suggested by 
iKaufi^mann elall (l2003l) forz < 0.1. 

4.2. High redshift diagnostic diagram (0.548 <z< 0.884J 

For the objects placed at higher redshifts, where we lack useful 
AGN-activity indicators such as the [Nil], [SII] or [OI] Unes, a 
different kind of diagnosti c diagram is needed to determine the 
galaxy spectral type. The iLamareille et al.l (l2004h involves the 
[Oil] line (Fig. |5j and represents a high-redshift option to the 
more classic red diagnostic diagrams. Unfo rtunately, this dia - 
gram is strongly biased against composites (lLamareillell2010h . 
which overlap with star-forming ga laxies and Low Ionization 
Nuclear Emission Regions (LINERs. lH"eckmanll980h . The [Oil] 
emission line is thought to be either an indicator of ongoing star 
formation or AGN-activity, due to the relatively low ionization 
potential that is required compared to the [Nil] or [SII] transi- 
tions. As a consequence, this diagnostic diagram cannot be con- 
sidered as reliable as the set of low-redshif t diagnostic diagrams 
(iBongiorno et al.ll2010l:lJuneau et al.ll201 ih . 
In this diagram it becomes necessary to correct line ratios for 
reddening. This is because the emission-lines that are used as 
indicators are placed far away from each other in the spectra 
and they are thus differently affected by reddening. Another 
way to overcome this problem is to use the EW of the emis- 
sion lines. The use of the EW instead than the extinction cor- 
rected fluxes does not change our galaxy classification signifi- 
cantly (see Appendix B). 

For separating AGN from star-forming galaxi es (Fig. |5}, we used 
lines which have been empirically defined bv ILamar eille (20 1 Cft 
by using galaxies from the SDSS. The left panel of Fig.[5]shows 
the distribution of SDSS galaxies (in gray) and the zCOSMOS 
galaxies superimposed on it. For every redshift bin (indicated 
with different colors) a higher mass means, in most of the cases, 
a shift towards higher values of [OII]/Hj6. Although the stacked 
galaxies show, as in the [Nil] diagram (Fig.lUi, a trend with the 
total stellar mass, almost all the objects are placed in the SF zone. 
For the higher masses, we notice that some of the stacked galax- 
ies fall in the region of mixed contribution from star-forming 
galaxies and LINERs, while only the points representing the 
highest mass-redshift bin falls in the star-forming/Seyfert region. 
This poor ability to properly classify AGN using the [OII]/H;6 
ratios was also noticed in the analysis of z COSMOS spectra for 
24 micron-selected galaxies dCaputi et al.l i2008. 2009). The di- 
agram is known to be the most sensitive to shocks. This might 
be the reason why some objects are likely to be misclassified as 
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Fig. 6. The new diagnostic diagram (MEx) by Juneau et al. 
(2011) makes use of the stellar mass instead of the [Nil]/!!!!- 
emission line ratio. In this plot, every color represents the evolu- 
tion in mass for a fixed redshift bin (all the z range): 0. 177 < z < 
0.242 (black), 0.242 < z < 0.306 (purple), 0.306 < z < 0.371 
(blue), 0.371 < z < 0.436 (light blue), 0.548 < z < 0.632 
(green), 0.632 < z < 0.716 (yellow), 0.716 < z < 0.8 (orange) 
and 0.8 < z < 0.884 (red). The mass increases from left to right 
along each track. The black diamonds indicate mass incomplete- 
ness (see Fig. [U. 



star-forming galaxies rather that LINERs, which would appear 
in a similar position along the y axis, but at higher [OII]/H;6 val- 



4.3. New generation of diagnostic diagrams 

Classical diagnostic diagrams are largely used by the scientific 
community to classify AGN from emission line ratios, especially 
in the optical range. Nevertheless, these diagrams present some 
problems that make it preferable to look for new and mo re effec- 
tive diagnostic methods. Recentlv. Ijuneau et alJ (1201 ih adopted 
a new diagram (MEx diagram) that can be used at higher redshift 
with respect to the classical diagnostic diagrams. A correlation 
between galaxy stellar mas s and metal content ha s been observed 
(mass-metallicity relation, iTremonti et al.ll2004t ISavaglio et al.l 
20(3) and the [NII]/Ha ratio is known to trace metallicity for 
SFGs dKewlev et all 120061: IStasiiiska et alJl2006h . For this rea- 
sons, the authors chose to substitute the emission line ratio 
[NII]/!!^ in the BPT diagram with the total stellar mass. This 
new kind of diagnostics has the important advantage of allowing 
to classify all the objects of a quite wide redshift range within 
a unique diagram. Moreover, it avoids some of the issues linked 
to the reliability of the [Oil] diagram, though the MEx diagram 
suffers from more blending between LINERs and Seyfert 2s with 
respect to the latter 

The mass trend we akeady observed in the [Nil] diagram is now 
clear in the MEx diagram (Fig.|6]l, where we can represent galax- 
ies spanning our entire redshift range. As in Fig. |4]and FigJSJ 
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each color represents a different redshift bin and tfie points indi- 
cate the mass bins. At log M»/Mq > 10.2, the galaxies leave the 
SF region to enter the composite region of the diagram. With the 
increase of the mass, the [OIII]/H;6 ratio decreases up to the point 
where the galaxies start to be classified as composites. Then, 
the ratio increases again, as already noted in the [Nil] diagram. 
While the stellar mass tracks overlap in Fig. |4]and FigjS] they 
are parallel to each other in the SF region of the MEx diagram 
and systematically offset towards the AGN region for increasing 
redshift (Fig. |6]l. Therefore, higher the redshift higher the num- 
ber of composites which are classified as transitional objects or 
AGN. In this respect, the MEx diagram is the one that shows the 
most clear sign of redshift evolution. 

5. Discussion 

In this section, we discuss our main results and explore future 
possibilities. 



5.1. Average and median stacking 

Stacking analyses may be challenging to interpret if the stacked 
galaxies span a broad range of physical properties. When taking 
averages, galaxies with strong lines might dominate the signal 
over objects with intrinsically weaker lines. As an example, we 
would expect massive galaxies that are star-forming and LINERs 
to have much weaker lines (in particular [OIII]) compared to 
Seyferts with the same stellar mass. Therefore, a potential risk 
would be that a few percent of Seyferts dominate the signal and 
raise [OIII] on the stacked spectrum, yielding a biased view that 
the stacked galaxies are AGN dominated by number In order 
to search for this effect, we compared the results with median 
stacking and average stacking. We made the comparison for a 
sub-sample of two median stacked redshift bins (15 stacks in to- 
tal). The line ratios we calculated from the stacks obtained with 
median and average stacking are very similar and lead to the 
same spectroscopic classification. Therefore, we did not find any 
indication of a bias towards more powerful AGN dominating the 
signal of the stacks. We conclude that there is no substantial dif- 
ference between median and average stacking of galaxies in our 
mass-redshift bins. All the measurements shown in this paper 
are from average stacked spectra. We can identify two possible 
reasons for the similarity between the measurements: a) From 
the definition of average and median stacking, we expect the dif- 
ference between the two to be noticeable only when the num- 
ber of AGN in a particular mass-z bin is very low (compared 
to the number of star-forming galaxies). From the results which 
are shown in the paper and as already expected from previous 
studies, this is more likely to happen at low masses. However, 
the chance of noticing a difference between median and average 
stacking is reduced in case of symmetric shape of the underlying 
distribution (e.g. Gaussian), when in a large sample the median 
value of the distribution tends to the average value. In our case, 
we deal with hundreds of galaxies per bin, thus with a statistics 
which might be large enough to hide this effect, b) The effect is 
still there, but it is not possible to notice the difference between 
median and average stacking and the trend is somehow lost in 
the measurements errors, which we estimate to be < 20% (see 
Section 3.4). The typical difference we find on the emission-line 
ratios is, in the [Nil] diagram, < 5% and - 20% for the [NII]/!!^ 
and [OIII]/Hy6 ratios, respectively. In the [Oil] diagram, the dif- 
ference is ~ 10% and ~ 40% for the [OII]/H/3 and [OIII]/H/3 
ratios, respectively. 
However, we should keep in mind that in our study all the low 
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Fig. 7. Measured E(B-V) (inferred from the Balmer ratios) as a 
function of stellar mass, with tracks corresponding to a given 
redshift: 0.177 < z < 0.242 (black), 0.242 < z < 0.306 (pur- 
ple), 0.306 < z < 0.371 (blue), 0.371 < z < 0.436 (light 
blue), 0.548 < z < 0.632 (green), 0.632 < z < 0.716 (yellow), 
0.716 < z < 0.8 (orange) and 0.8 < z < 0.884 (red) (see Fig.©. 



redshift mass-redshift bins and several of the high redshift bins 
are considered to be complete (see Sect. 2.1, Sect. 3.1). Hence, 
each stack is representative of the average behavior of a particu- 
lar bin. 



5.2. Visual extinction 

It is known that dust extinction in star-forming galaxies depends 
upon star f ormation rate (SFR) , metallicity and stellar mass inde- 
pendently (iGarn & Bestll2010b . High extinction values are asso- 
ciated with an increase of any of these three parameters, although 
the fundamental property connected with extinction seems to be 
the stellar mass. However, in our study we consider galaxies of 
different spectral type. We would expect to find increasing values 
of the visual extinction for higher masses, since the high mass 
tail of the galaxy mass function of galaxies placed at high red- 
shift (z ~ 2) is populated by Ultra Luminous Infrared Galaxies 
(ULIRG s, infrared lumino sity at 8 - 1000 microns higher than 
10'^ Ln. lDaddi et alj|2005[) that store a large amount of dust red- 
dening the outcoming light. On the other hand, we may expect 
the opposite trend, namely a decreasing value of the visual ex- 
tinction with the stellar mass, since on the high tail of the galaxy 
mass distribution we observe objects that might have experi- 
enced the AGN-feedback (e.g. jets blowing up the gas and dust 
content). In general, massive systems that have already experi- 
enced star formation in the past are found not to contain a large 
amount of gas and dust. Less massive galaxies, according to the 
downsizing scenario, are more likely in the evolutionary stage 
where they still have material to fuel the star-formation. 
We find that the distribution of the E(B-V) values is nearly flat 
for increasing mass and redshift. In particular, the visual ex- 
tinction measured from the stellar-continuum subtracted spectra 
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Fig. 8. Classical [SII] diagnostic diagram. In both panels demarcation curves bv lKewlev et alj ( l2006h separate star- forming galaxies 
from Seyferts and LINERs. Each color represents the trend with stellar mass for a fixed redshift bin: 0.177 < z < 0.242 (blue), 
0.242 < z < 0.306 (green), 0.306 < z < 0.371 (purple) and 0.371 < z < 0.436 (red). T he total stellar mass increases along each 
track, starting from the top. In the left panel, the SDSS galaxies from lVitale et al.l (1201 2|) are plotted in gray. The right panel offers 
a closer view of the zCOSMOS galaxies. 



shows an increment with the stellar mass only in the low -redshift 
bins (Fig.|7]and Tab.|2]i. At higher redshifts the same trend may 
be hidden by the bigger errors in the measurement of E(B-V) 
(see Sect. 3.4). On the other hand, the E(B-V) values inferred 
from the SED fitting results displayed in Tab. |2]show a different 
tentative trend. The E(B-V) values are higher in the high-redshift 
bins, but with no distinction between high and low masses. These 
evidence are not sufficient to indicate what the main contribution 
to extinction could be - reddening from dusty ULIRGs or low 
extinction due to AGN-feedback. The results seems to weakly 
suggest a higher extinction for the more massive galaxies, espe- 
cially if placed at higher redshift. However, further studies are 
needed to shed light on this topic. 

5.3. Diagnostic diagrams 

Emission lines are a good indicator of star- forming or AGN ac- 
tivity. By means of diagnostic diagrams and spectral stacking, 
we found the indication for a link between the total stellar mass 
and the chance of identifying galaxies as AGN. A temporal evo- 
lution - namely, an observed trend with the redshift - is also pos- 
sible but, even if present, it is not well visible within the classical 
diagnostic diagrams. This could be due to our redshift bins being 
too close to each other for the stacked galaxies to show a trend 
with redshift in the classification provided by the diagrams. It is 
important to notice that we are studying a cosmic period (z < 1) 
that is considered as crucial for AGN-feedback and quenching of 
the star formation to happen, as it has seen both the drop in the 
star form ation history and the build up of a bout half of the red 
sequence dBell et alj2007tlFaber et al.ll2007h . However, the time 
interval over which to probe evolutionary effects is larger when 
using the MEx diagram (0.17 < z < 0.88) than when using the 
BPT (0.17 < z < 0.44) or [Oil] (0.55 < z < 0.88) diagrams 
alone. 
The indication of an evolutionary trend is present in the MEx 



diagram (Fig. |6]l. Here the data points representing galaxies 
in the same mass range are progressively offset to higher val- 
ues of [OIII]/Hj6 for increasing redshift. On the other hand, 
along each track (fixed redshift, increasing stellar mass) the data 
show a metal enhancement (see also Fig. |4|i. The reason for 
the latter could be researched in the mass-metallicity relation 
and the efficiency of galactic winds in r emoving metals from 
low-mass galaxies ('L equeux et al.l fl979t iTremonti et all 120041 : 
Maiolino et al. 200 8). Outflows generated by starburst winds 
may easily eject metal-enriched gas from low-mass galaxies 
placed at low redshifts, m aking their enrichmen t less significant 
than in massive systems ( ITremonti et al.ll2004b . Moreover, low 
mass systems are, ac c ording to the galaxy downsizing scenario 
(iJuneau et alJ l2005t lAsari et al.l 120071: iPerez-Gonzalez et al.l 
|2008|) at an earlier evolutionary stage, where they are still con- 
verting most of their gas into stars. For this reason they are 
poorly metal-enriched compared to massive galaxies. However, 
it should be kept in mind that our purely magnitude-limited sam- 
ple suffers from selection biases and mass incompleteness to- 
wards the highest redshifts, where most of the AGN are found 
among the stacked galaxies. The set of classical diagnostic di- 
agrams turned out to be not completely suitable for our study. 
While the [Nil] diagram nicely shows that the mass can deter- 
mine the galaxy classification, the [SII] - due to the problems al- 
ready discussed in Sect. 4.1- does not show the same trend and 
classifies all the stacked galaxies as SF (Fig.jS]). The [OII]-based 
diagram, on the other hand, leaves us with big uncertainties be- 
cause of the intrinsic nature of the diagram itself (rather s ensitive 
to shocks, it is not able to disentangle the contribution of stellar 
sources from AGN). As a matter of fact, the diagram shows an 
overlap between star-forming galaxies and Seyfert 2s, as well 
as an overlap between the former and composites. However, the 
higher-mass galaxies are closer to the demarcation line and the 
AGN region. The usage of the MEx diagnostic diagram partially 
solves the issues related to the ineffectiveness of the classical 
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diagnostic diagram at high redshift ([Oil]) and serves to the pur- 
pose of looking for AGN contamination at high-masses (see next 
section). The mass threshold above which the objects start to be 
classified as Composite and AGN is well visible. 

5.4. Quenching of the star-formation 

We noticed a general increase of the age of the stellar popula- 
tions for increasing mass, at both high and low redshifts. Tab. 
|2] summarizes all the information on the stacked galaxy spec- 
tra in the stellar mass and redshift bins illustrated in Fig. [1] The 
most suitable stellar continuum template to be subtracted to each 
stacked spectrum has been chosen from the set of 234 stellar 
templates (39 BC templates x 5 reddened spectra). The table 
shows the reddening inferred from the stellar template, E(B-V)j, 
and the reddening calculated from the Balmer lines, E(B-V)„,, 
whose fluxes have been measured from the stellar-continuum 
subtracted spectra. The metallicity of the stellar template gives 
an import ant indication o f the presence of the age-metallicity de- 
generacy dWorthevI 19941) . For a fixed redshift and along the mass 
sequence (Fig. [1] Tab. |2]), the stellar populations become pro- 
gressively older A young stellar population following an older 
one in the mass sequence presents, in most of the cases, a higher 
metallicity and a higher reddening as inferred from the stellar 
fit. A stellar population following one of the same age along the 
mass sequence is often characterized by equal or higher values 
ofZandE(B-V),. 

This evidence finds a possible explanation in the mass-assembly 
downsizing scenario, where most of the mas sive galaxies build 



their mass earher than lower mass galaxie s (jCo wie et al.l 



Kodama et al] 12004 iThomas et al.1 120051: Ffanakaet all 12005 : 
Juneau et al.l l2005t iBundv et al.ll2006h . According to this sce- 



is galaxie s (ICo \ 
aiTl2005l: Ffana 



1996; 



nario, the number of red massive objects is constant up to z = 1, 
dominating as the main contribution to the mass population, 
while the number density of red galaxies with M < 10" M© 
increases with cosmic time. Conversely, the number density 
of blue galaxies decreases from z = 1 to today. Downsizing 
suggests that something turns the galaxy transformation out at 
around z - 1 . If the most massive objects build their mass dur- 
ing an earlier period of time with respect to the less massive 
galaxies, we expect to find older stellar populations in the for- 
mer because they had more time to evolve. 
If the effect of having younger stellar populations at higher red- 
shifts was already expected from different galaxy evolutionary 
scenarios, the evidence of older stellar populations with increas- 
ing mass and fixed redshift points to the downsizing scenario in 
particular It is likely that, in these objects, a mechanism acted 
to quench the star formation. Now the question is: can this pro- 
cess be an AGN-feedback? Considering the downsizing scenario 
where young low-mass galaxies are still undergoing star forma- 
tion and high mass galaxies are the product of a earlier mass 
assembly, the latter are passively-evolving galaxies with no star- 
formation episodes. In this respect, AGN might have just sup- 
pressed star formation with their feedback, by means of trans- 
forming blue galaxies into red ones that continue to accrete 
mass. As a support to this theory, it has been found that the stel- 
lar populatio n of the host gala xies appears older after the AGN 
active phase dSarzi et al.ll2005h . Nowadays, several theories and 
models include the AGN-feedback, but observational evidence 
are still needed to prove them right. 
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Table 2. Information on the stacked galaxies of the mass-redshift bin of Fig. [T] 







SSP 2.5 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


SSP 2.5 Gyr 


0.2 


SSP 2.5 Gyr 


0.2 


SSP 1.4 Gyr 





SSP 1.4 Gyr 


0.2 






AGN 




AGN 




AGN 




AGN 




AGN 




AGN 




AGN 








0.451 


0.08 


0.415 


0.5 


0.011 


0.5 


0.153 


0.08 


0.333 


0.08 


0.369 


0.5 


0.459 


0.2 


SSP1.4Gyr 


0.2 


SSP 2.5 Gyr 


0.2 


SSP 2.5 Gyr 


0.2 


SSP 11 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


AGN 




AGN 




AGN 




Comp 




AGN 




AGN 




AGN 




AGN 




0.613 


0.5 


0.437 


0.2 


0.580 


0.08 


0.687 


0.2 


0.317 


0.2 


0.240 


0.2 


0.023 


0.08 


0.125 


0.08 


SSP 2.5 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


SSP 2.5 Gyr 


0.2 


SSP 900 Myr 


0.2 


SSP 1.4 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


Exp 12 Gyr 


0.2 


SSP 900 Myr 





Comp 




Comp 




Comp 




Comp 




Comp 




Comp 




AGN 




AGN 




0.539 


0.08 


0.363 


0.2 


0.580 


0.08 


0.240 


0.5 


0.196 


0.08 


0.449 


0.08 


0.331 


0.08 


0.170 


0.5 


Exp 12 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


SSP 900 Myr 


0.2 


SSP 900 Myr 


0.2 


CST 6 Gyr 


0.4 


CST 6 Gyr 


0.4 


Exp 12 Gyr 


0.2 


Exp 12 Gyr 


0.2 


SF 




Comp 




Vomp 




SF 




Comp 




Comp 




Comp 




Comp 




0.389 


0.2 


0.403 


0.2 


0.318 


0.5 


0.342 


0.5 


0.635 


0.08 


0.349 


0.08 


0.074 


0.08 


0.235 


0.08 


SSP 900 Myr 


0.2 


SSP 900 Myr 


0.2 


SSP 1.4 Gyr 


0.2 


SSP 1.4 Gyr 


0.2 


CST 6 Gyr 


0.4 


CST 6 Gyr 


0.4 


SSP 25 Myr 


0.6 


CST 6 Gyr 


0.2 


SF 




SF 




SF 




SF 




Comp 




Comp 




Comp 




Comp 




0.365 


0.5 


0.369 


0.5 


0.465 


0.08 


0.391 


0.08 


0.505 


0.08 


0.350 


0.08 


0.563 


0.08 


0.443 


0.2 


Exp 12 Gyr 


0.2 


Exp 12 Gyr 


0.2 


SSP 640 Myr 


0.2 


SSP 1.4 Gyr 





CST 6 Gyr 


0.2 


CST 6 Gyr 


0.2 


CST 6 Gyr 


0.2 


SSP 25 Myr 


0.4 


SF 




SF 




SF 




SF 




SF 




SF 




SF 




^.omp 




0.340 


0.2 


0.384 


0.2 


0.322 


0.5 


0.244 


0.2 


0.341 


0.08 


0.373 


0.08 


0.567 


0.2 


0.401 


0.08 


SSP 640 Myr 


0.2 


SSP 640 Myr 


0.2 


SSP 900 Myr 


0.2 


Exp 12 Gyr 





CST 6 Gyr 


0.2 


SSP 25 Myr 


0.4 


SSP 25 Myr 


0.4 


CST 6 Gyr 





SF 




SF 




SF 




SF 




SF 




SF 




SF 




SF 




0.264 


0.5 


0.298 


0.5 


0.231 


0.08 


0.224 


0.08 


0.362 


0.08 


0.172 


0.08 


0.357 


0.08 


0.197 


0.5 


SSP 290 Myr 


0.4 


SSP 640 Myr 


0.2 


Exp 12 Gyr 





SSP 640 Myr 
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SF 




SF 




SF 




SF 




















0.292 


0.2 


0.188 


0.2 


0.253 


0.08 


0.182 


0.08 


















SSP 900 Myr 


0.2 


SSP 1.4 Gyr 













SP 


E(B-V), 


SF 




SF 
























Class 




0.167 


0.08 


0.345 


0.08 






















E(B-V)„, 


Z 



Notes. In each box: template of stellar population that fits the stacked zCOSMOS spectrum the best (SP, top-left); reddening of the stellar spectrum (E(B-V)j top-right); value of the extinction as 
measured from the Balmer ratios after stellar continuum subtraction (E(B-V),„, bottom-left); metallicity (Z, bottom-right); spectral classification (Class, middle of the box) according to the MEx 
diagnostic diagram. SSP, CSP and Exp stand for Simple Stellar Population, Constant Stellar Template and Exponential decay of the star-formation, respectively. The legend is to the bottom-right of 
the table. 
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5.5. AGN identification 



The downsizing scenario takes into account a quenching mech- 
anism to explain galaxy evolution. Our findings fit into the big 
frame of this latter scenario. Tab. |2]reports the spectral classifica- 
tion (AGN, Composite or SF) of the stacked zCOSMOS galaxies 
according to the MEx diagnostic diagram (Fig. |6). It is clearly 
visible that AGNs are common amongst the high-mass bins. This 
is particularly the case of higher redshift bins. However, some 
stacked spectra are placed right on the demarcation curves that 
separate AGN from Composites, or right beside them. In these 
cases, it might be either that the stacking galaxies do not belong 
to a single population - namely, they are a mixture of AGN and 
star-forming galaxies - or that the underlying population is truly 
dominated by composite galaxies. 

The use of stacked spectra instead of individual spectra leads to 
some difficulties in estimating the number of SFGs and AGN. If, 
from one hand, this technique allows to study the average proper- 
ties of galaxies in a given mass-redshift bin - providing average 
spectra for otherwise individually unclassifiable galaxies with 
too low S/N - from the other hand it does not provide us with the 
exact number of spectral types. Furthermore, each bin collects 
a variable number of galaxies, due to the presence of filaments 
(overdensities), voids (underdensities) and selection effects act- 
ing at the highest as well as at the lowest stellar masses (Fig. [T| 
and Tab. [TJ. Therefore, it is not possible to compare the exact 
number of AGN between the different bins. Assuming that all 
the galaxies in a given stack have the same classification as the 
stacked spectrum, we obtain the following statistics. The number 
of single galaxies that are spectroscopically identified as AGN in 
the high-redshift bins (5 568 galaxies in total) is 1 690 (30.3%), 
while the Composites are 2 491 (44.7%). This statistics is higher 
compared to the galaxies placed in bins at lower redshift (3 742 
galaxies), where we find 395 AGN (10.5%) and 775 Composites 
(20.7%). However, assuming that all the galaxies in a given com- 
posite or AGN bin are indeed AGNs will likely result in an over- 
estimate of the true AGN fraction. Nevertheless, the relative dif- 
ference between the higher and lower redshift bins remain an in- 
teresting hint of a greater AGN fractio n at earlier epoch. B ased 
on the study of individual galaxy spectra llvezic et al1(l2002h find, 
with SDSS data, a starburst over AGN ratio equal to 18, which 
is a much higher ratio than what we find. 

The search for AGN could benefit from the complementary 
analysis of individual galaxy spectra or from analyses at other 
wavelengths. The X-domain can identify AGN because the lu- 
minous, compact X-ray emission {L2-10 kev > 10'*^ ergs s"') 
is an almost certain indicator of the existence of an AGN 
(iBauer et al.ll2004t iBrandt & HasingeJl2005l: IComastri & Brusal 



l2008t iBrandt & AlexandeJ 1201 Oh . given the extremely small 
contribution from star formation to the overall emission at these 
frequencies. X-ray observations have been found to be very ef- 
ficient in revealing accreting black holes in galaxies which were 
not classified as AGN from the analysis of optical data (see 
iBrusaet aDl2010l and references therein). This suggests that 
the use of multi-wavelength data can be beneficial in AGN-host 
studie s. The XMM-Newton wide-field survey in the COSMOS 
field (iHasinger et al.l l2007h provides a large sample of point- 
like X-ray sources (~ 1 800) with complete ultraviolet to mid- 
infrared (including Spitzer data) and radio coverage. The sur- 
vey - thanks to the high efficiency of X-ray observations in 
identifying AGN - currently allows further investigations (e.g. 
iBrusa et al]|2010l) and contributes to answering questions on 
galaxy-black hole coevolution up to high redshifts. 



5.6. Aperture effects 

Aperture biases are particularly important in emission-line stud- 
ies. Spectroscopic observations are performed using different 
aperture sizes, depending on the specific instrument. This is- 
sue plays a critical role in disentangling the main contribution to 
galaxy emission. The so-called "aperture effect" tells us that the 
shape and size of the instrument aperture influences the galaxy 
classification. We would expect big fibers/slits to select emis- 
sion from all over the galaxy, hence several contributions to the 
line emission (AGN from the nucleus, starburst from the outer 
regions of the host galaxy) are present. On the other hand, a 
smaller aperture would mainly select light coming from the nu- 
cleus and the bulge, missing the disk where most of the star for- 
mation takes place. For example, SDSS data are known to be 
affected by aperture effect due to the large fixed size (3") of the 
fiber used for observations which includes on average > 20% 
of galaxies area at z> 0.04 (iKewlev et alj|200q) . Therefore, the 
stellar contribution to the emission from the host causes a large 
number of sources to be classified as transitional objects in the 
diagnostic diagrams. The same argument applies at higher red- 
shift, where the slit width becomes comparable to the typical 
angular size of galaxies. Furthermore, higher-redshift galaxies 
present higher star-formation rates, so the effect of host galaxy 
dilution of the AGN signatures could be even more important 
and yield more often composite signatures. From this point of 
view, with the 1" width of the slits used for the zCOSMOS 
observations and a seeing always better than 1.2" (iLillv et al.l 
l2007h . our sample is less sensitive to contamination due to emis- 
sion coming from the most outer regions of the galaxies with 
respect to other studies at similar redshift. 
However, the conventional belief that smaller apertures select 
more AGN-like emission has been challenged by Shiel ds et al.l 
(l2007h . The authors show that for many co mposite objects id en- 
tified in the Palomar spectroscopic survey (iHo etal.ll 1 997bl) the 
observed line-ratios do not appear more AGN-like going on 
smaller apertures (10 - 20 pc). Therefore, the aperture effect is- 
sue remains controversial and it is worth further studies. 



6. Summary and conclusions 

In this work, we used the zCOSMOS 20k catalogue to create a 
sample of galaxies with reliable spectroscopic redshift and with- 
out broad-line AGN, in order avoid overshining issues (AGN 
brighter than the surrounding galaxy) and investigate the evolu- 
tion of galaxy spectral properties up to z ~ 1 . The sample, con- 
tainig ~ 1.5 x 10'* objects, has been divided in 61 bins over the 
stellar mass-redshift plane (Fig.[TJ. This was done taking into ac- 
count the completeness in mass and the gap in redshift, this latter 
due to the impossibility of measuring all the emission lines of the 
classical diagnostic diagrams. In order to improve the signal to 
noise (S/N) of the single spectra, we have combined all the spec- 
tra in each bin to obtain stacks, allowing accurate flux measure- 
ments. After stacking the spectra, we performed the stellar con- 
tinuum subtraction on the galaxy composites, using templates of 
stellar population synthesis. This provides new spectra that are 
almost free from absorption contaminations by stars, and also 
gives information on the age of the stellar populations and their 
composition (metallicity). Finally, we applied optical emission- 
lines diagnostics to search for AGN contamination at the highest 
stellar masses. 
Our main findings are: 

- The visual extinction of the galaxies inferred from the 
Balmer ratios does not show a strong trend with the total 
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Stellar mass all over the z range. It increases slightly at low 
redshifts for increasing mass, but there is no indication of a 
clear trend at higher redshifts. 

- Galaxy stellar populations are older (in terms of time from 
the last burst of star formation) for more massive objects. 
This is in agreement with an evolutionary scenario that ac- 
counts for a fast mass assembling and stellar population ag- 
ing in case of massive objects. 

- In general, galaxy stacked spectra are more likely to be clas- 
sified as AGN rather that star- forming galaxies for large val- 
ues of their total stellar mass. In particular, we found that 
galaxies with log M,,/A^0 > 10.2 start to be classified as com- 
posite objects (where both photoionization by stars and nu- 
clear activity contribute to produce the emission lines). They 
are placed in the AGN region of the diagnostic diagrams for 
the highest masses considered in our sample. 

- At fixed redshift, the classification of the stacked spectra dis- 
plays a trend with increasing stellar mass. The general trend 
follows the left branch of the SDSS classification down to the 
bulk of the star-forming galaxies population. At the highest 
stellar masses, the tracks fall in the composite and then in the 
AGN region. This is visible in the [Nil] and in the MEx dia- 
gram, whereas it is not well observed in the [Oil] diagram. 

- The MEx diagram represents an useful tool to investigate the 
ionizing mechanisms inside galaxies, especially at high red- 
shift. The Lamareille (2010) [Oil] diagram is not as effective 
due to the ambiguity between star-forming, composites and 
Seyfert 2 galaxies in some regions of the diagram. Moreover, 
there is the need to be in a specific redshift interval to mea- 
sure all the involved emission lines. The confusion that this 
diagram shows can be due to its better sensitivity to shocks 
than to AGN photoionization. 

- While there is no clear trend with redshift on the low -redshift 
BPT diagram, there is a trend on the MEx diagram when 
combining the low- and high-redshift bins. It is visible that 
the high-redshift stacks display more composite- and AGN- 
like spectra. However, the explanation of this trend (higher 
AGN-detection rate at higher redshift) may include true evo- 
lutionary effects as well as selection biases. 

- The link between stellar population age and galaxy stellar 
mass, combined with the increasing AGN-detection rate for 
increasing mass and redshift, is consistent with a scenario 
where AGN could act to quench the star formation and, then, 
contribute to the transformation from young blue late-type 
galaxies to old red early-type galaxies. 
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Appendix A: 



Here we present the zCOSMOS stacked spectra which have been 
used in our study. 

For a fixed redshift, the stacked spectra - before stellar con- 
tinuum subtraction - show a [OIII] line whose strength decreases 
with the progressive increase of the stellar mass (Fig. lA.l) . es- 
pecially if compared with the H/3 line. The latter shows a weak- 
ening as well. The Ha line is very strong at low masses. The 
intensity of the [Nil] and Ha decreases up to the point, at the 
highest mass bin, where these lines show the same intensity. In 
Fig |A.2l the [OIII] and Hfi lines show the same trend as at lower 
redshifts. The [Oil] line decreases with the increase of the mass, 
but it is always the strongest line in the high-redshift composite 
spectra. 



Appendix B: 



In this work we used the dust-corrected flux ratios for [OII]/H 6 
instead of the EW ratios originally used bv iLamareillel (1201 Oh . 
This appendix (Fig lB.lb shows that the dust-corrected fluxes cor- 
respond well enou gh to the EWs for a sample of SDSS galaxies 
(I Vitale et alJl2012|) . especially in the regions occupied by our 
sam ple. Hence we can use the original demarcation lines from 
iLam areille (2010). However, the dust-corrected flux ratios span 
a slightly more restricted range of [OIII]/Hy6 values than the EW 
ratios. 
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